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Abstract. 

We revisit the study of the emptiness formation probability, the probability of forming a sequence of 
£ spins with the same ferromagnetic orientation in the ground-state of a quantum spin chain. We focus 
on two different examples, exhibiting strikingly different behavior: the XXZ and Ising chains. One has a 
conserved number of particles, the other does not. In the latter we show that the sequence of fixed spins 
can be viewed as an additional boundary in imaginary time. We then use conformal field theory (CFT) 
techniques to derive all universal terms in its scaling, and provide checks in free fermionic systems. These 
are based on numerical simulations or, when possible, mathematical results on the asymptotic behavior of 
Toeplitz and Toeplitz+Hankcl determinants. A perturbed CFT analysis uncovers an interesting l~ l \ogl 
correction, that also appears in the closely related spin full counting statistics. The XXZ case turns out 
to be more challenging. We use a simple qualitative picture in which the ferromagnetic sequence of spins 
freezes all degrees of freedom inside of a certain "arctic" region, that we determine numerically. We also 
provide numerical evidence for the existence of universal logarithmic terms, generated by the massless field 
theory living outside of the arctic region. 
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1. Introduction 

1.1. Emptiness formation probability in spin chains 

In recent years, a great deal of effort has been put into the understanding and calculation of correlations in 
integrable one-dimensional quantum critical systems. Although many physical properties can be understood 
from the Bethe-ansatz solution [1, 2], the explicit calculation of physical correlation functions is a much more 
formidable task. Much progress has however occurred in this direction, in particular in the spin-1/2 XXZ 
spin chain, a key physical example where such methods are successful. New methods have been developed 
to compute such correlations and their asymptotics exactly, starting from the lattice model [3, 4, 5, 6, 7]. 
Moreover, the use of conformal field theory (CFT) [8] provides a way to access various universal exponents 
in a simple and transparent way. However, it does not give access to non-universal quantities, and usually 
still requires some input from more exact methods [9] . These two complementary approaches have also been 
shown to agree in several studies, see e.g. [10, 11]. 

In the first framework, arguably the simplest quantity one can compute is the probability of forming a 
sequence of £ consecutive spins with the same up orientation, the so-called emptiness formation probability 
(EFP). For this reason it has already been subject to quite extensive studies [12, 13, 14, 15, 16, 17, 18, 19, 20], 
especially (but not always) in the simpler limit of an infinite chain. In particular, many exact results have 
been derived, expressing the EFP using multiple integral or determinantal representations, from which an 
asymptotic expansion can sometimes be extracted. Note that the study of the EFP is also interesting in - 
simpler - free fermionic systems. Indeed in this case the EFP usually boils down to a Toeplitz determinant, 
whose interesting mathematical properties have attracted much attention since the celebrated Kaufman- 
Onsagcr spontaneous magnetization problem [21, 22]. 

Two intriguing features stand out in the various asymptotic results obtained for the EFP. First, the 
leading term can depend on the model. In logarithmic form the EFP is proportional to I 2 for the XXZ chain 
[14, 15], while it only scales proportionally to £ for the Ising chain [23]. This behavior can be traced back 
to the fact that the XXZ chain has a £7(1) symmetry, that is broken for Ising: in fermionic language the 
particle number is conserved in one case, but not in the other. Second, the EFP is sensitive to criticality. 
When the system is critical, power law corrections - logarithmic terms for the logarithmic EFP - to the 
leading term start to appear. While these terms can often be explained at the purely mathematical level (e.g. 
singularities in the symbol generating the Toeplitz determinant), they have so far remained unexplained at 
the field-theoretical level. 

This paper is a partial attempt at filling this gap. We show that the Ising case can be understood in 
full detail using CFT. As a consequence, we confirm the universality of the power law correction derived 
by Franchini and Abanov [23], and show that it is a simple fraction of the central charge of the CFT. We 
provide several generalizations to open systems and finite geometries, where we compute a universal scaling 
function, in the spirit of entanglement entropy studies. All these results are extensively checked in free- 
fermionic systems, combining numerical evaluations and asymptotic results on Toeplitz and Toeplitz+Hankel 
determinants. Following [25], we then carry out a perturbed CFT analysis of the subleading corrections, 
and show how the first goes slightly beyond the known Toeplitz results. As we shall see it takes a form 
proportional to £~ x log£, with a coefficient determined by the central charge, the geometry of the system, 
and a single non universal short distance cutoff, the so-called extrapolation length [24]. We also study 
several related quantities, including the Renyi-Shannon entropy and full counting statistics, where this type 
of subleading correction makes an appearance as well. 

In contrast, the XXZ case turns out to be much more subtle. We are nevertheless able to gather some 
evidence for the universality of the power law corrections, and show how the general scaling may be explained 
in terms of an "arctic phenomenon" , in which all degrees of freedom are frozen in a region of area proportional 
to I 2 in imaginary time. Such behavior is familiar in the study of domino tilings on the Aztec diamond. In 
this picture a massless field theory, whose exact properties are unknown, lives outside of a "frozen region" 
that we determine numerically in various free fermionic examples. 
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1.2. Problem studied 



Let us consider the following quantum Hamiltonian, so called XY chain in transverse field. 

L 



i 



1-1/ 

2 



/«7- 



(1) 



We will study this Hamiltonian with periodic boundary conditions (c£ +1 — a", a 



x,y), as well as open 

boundary conditions (c£ +1 = 0, a = x,y). The system is in the ground-state and the spins are 
measured in the basis of the eigenstates of the crj . What makes this model particularly simple is that using 
a Jordan- Wigner transformation, it can be rewritten in terms of free fermions [32] : 

L 

II = > (c]c i+1 + v 4c\ +1 + h.cj -hJ2 (2c]c, - l) , (2) 
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where c 



t 

L+l 



for periodic boundary conditions, and c 



.t 

i+i 



for open boundary conditions. Here 



N = Y2j=i c j c j * s ^ ne fermion number operator. The dictionary between spins and fermions is fairly simple: 
a fermion at site j corresponds to an up spin, while a vacant site j corresponds to a down spin. 

As is well known, this model has several critical lines, as a function of the parameters v and h. The first 
we consider is v > and h = 1, and belongs to the Ising universality class]:. Its low-energy properties are 
described by the Majorana fermion, a CFT with central charge c = 1/2. In particular, the point h = 1, v = 1 
is the Ising chain in transverse field (ICTF). The other (y — 0, \h\ < 1) has a vanishing fermion pairing term, 
and therefore an additional U(l) symmetry. It belongs to the free compactified boson universality class (or 
Luttinger liquid), and has central charge c = 1. In the following, we will mostly focus on these two critical 
lines. 
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(a) Periodic chain 

Figure 1: Emptiness formation probability for a sequence of £ = 
L. For the open chain the fixed spins are the £ first. 



i L-l 
(b) Open chain 

8 consecutive up spins in a system of length 



One may then ask what is the probability of observing a sequence of £ consecutive up spins, somewhere in 
the chain (in the open chain, the sequence ends on the left edge). More precisely, let us denote by A the 
subsystem formed by the £ spins, and B the rest. See Fig. 1 for a picture. The reduced density matrix for 
A is 

p^TrshWI, (3) 
and the EFP is defined as 

7>Htt...t|/M|tt...t>. (4) 
For later convenience, we also introduce the logarithmic emptiness formation probability (LEFP) 

£ = -\ogV. (5) 

We wish to study the asymptotic behavior of £ with the size £, while keeping the aspect ratio fixed to a 
certain fixed value £/L G [0; 1). We are mainly interested in identifying universal terms in the process. In 
the following we use the notation £ p in the periodic system, and £ in the open system. 

X We ignore here the other similar line v > 0, h = —1, because it involves changes in conformal boundary conditions. We 
discuss this in Sec. 5.1. 
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1.3. Outline and summary of the results 

The paper is organized as follows. In Sec. 2 we study the example of the Ising critical line, and show that the 
ferromagnetic string acts as an additional boundary in euclidean time. We then use boundary CFT methods 
to compute universal terms in the LEFP. In particular we obtain 



£ v = a x t + - log 



£ Q = ail - — log 
16 



L . 




— sin 


;?)] 


7T 





4£ tan 2 (f^ 
7T sin (^) 



& „ 2-cos(f)logf (1 



X 



32 sin(f) L 



Oil). (7) 



c is the (universal) central charge, a%, Og , are non universal, we compute them using asymptotic results 
on Toeplitz and Toeplitz + Hankcl determinants in Sec. 3. The "extrapolation length" £ is also a non 
universal quantity, known in the context of surface critical phenomena [24]. We show that 

* = k (8) 

for the Ising critical line of the XY chain. The L^ 1 log£ terms can be seen as semi-universal [25]: once £ is 
fixed they are solely determined by the central charge of the CFT, and the geometry of the system. For the 
EFP this type of correction goes slightly beyond the available mathematical results on Toeplitz determinants, 
and offers an interesting generalization to a recently raised conjecture [26]. 

We then turn our attention (Sec. 4) to the XX(Z) chain, as an example of a free c = 1 boson CFT. This 
is known to exhibit strikingly different behavior, with in particular the LEFP scaling as 

£ P ,o = a 2 e 2 + a 1 £ + b \og£ + O(l). (9) 

We propose a simple qualitative picture for such a behavior, in which the ferromagnetic string in imaginary 
time generates a region of area oc £ 2 where the degrees of freedom are frozen. A massless held theory which 
is not strictly speaking a CFT, then lives at the exterior of this region. This picture is a manifestation of 
the well-known arctic phenomenon for dimers. We also provide numerical evidence for the universality of 
bo, by studying different lattice realizations and varying the aspect ratio £/L. The analytical derivation of 
bo using field-theoretical methods is however left as an important open problem. 

Finally in Sec. 5, we briefly discuss some related problems to which our results could be applied, including 
the Shannon entropy and full counting statistics. 

2. The ferroelectric string as a conformal boundary condition 

This section mainly relies on a simple observation: the configuration with £ up spins, viewed in euclidean 
time, should renormalize to a conformal boundary condition. There are only two possible (free or fixed) 
conformally-invariant boundary conditions for the Ising CFT [27]. The appropriate one can be determined 
by going back to the 2d classical Ising model. In a transfer matrix picture, the classical spins correspond, 
after taking the Hamiltonian limit, to the eigenstates of the a*. Therefore, a configuration with all fixed up 
spins in the a z basis translates into 

Itt ... t), = i (1-0* + !<->*) (I -Ox + !<->») ■ ■ ■ (1-0* + !<->*) (io) 



— T 



(11) 



= |free>* (12) 

a state where all orientations are allowed, with equal amplitude, in the x direction. Hence the configuration 
with up spins in the z basis should renormalize to a free conformal invariant boundary condition. 
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2.1. Imaginary time picture 



Let us now assume that the sequence of up spins can be considered as a free boundary condition in euclidean 
time. Since the external boundary conditions are also free, there are no changes in boundary conditions. 
The EFP is then given by the ratio of partition functions 

,(slit) 



P 

i) 



""cyl 

Z C y\ 



(13) 



Z cy i is the partition function of an infinite cylinder, and is the partition function of an infinite cylinder 

with a slit inserted in the middle, as is shown in Fig. 2a. In the open geometry the cylinders are replaced 
by strips, see Fig. 2b. In both cases, the logarithmic EFP is given by the corresponding differences in free 
energies F = — log Z . 



L 




(a) Infinite cylinder with a slit 



(b) Infinite strip with a slit 



Figure 2: Space-time geometry (the time direction is vertical). In the cylinder geometry, there are two 
corners with angle 2~k. In the strip geometry only one corner has angle 27r, but there are two additional 
corners with angle ir/2. 

Then, from general arguments, we expect the LEFP to scale as 

£ = - log 7? = a x l + b log L + f(£/L) + a Q + o(l) (14) 

Notice that the bulk free energies are canceled in the ratio (13). a± is the line free energy contribution of 
the slit. Let us now explain the next two terms. As was first pointed out by Cardy and Peschel [28], the 
logarithmic divergence is a general consequence of the presence of corners in the partition functions we have 
to evaluate. Indeed, there are two corners with angle 27r in , and one 2ix as well as two tt/2 corners in 



(slit) 



. <; ( ) . The coefficient b takes a simple form: for each corner with angle 9 there is a contribution 

r I H IT \ 7T 

log I 



AF = 



c 

24 



^bec 



(15) 



to the free energy, c is the central charge of the corresponding conformal field theory, with c = 1/2 for the 
Ising universality class. h{, cc is the dimension of a possible boundary changing operator, but there are none 
here (/i& cc = 0) since all boundaries are free. We therefore have bo = c/8 = 1/16 in the periodic geometry, 
and bo = — c/16 = —1/32 in the open geometry. f(£/L) is a universal function of the aspect ratio i/L, 
that we calculate in the next section, ao is a constant. Notice that both a\ and ao are non- universal and 
therefore not accessible within CFT. We calculate them for the Ising chain in Sec. 3, using exact results on 
the asymptotics of Toeplitz and Toeplitz + Hankel determinants. 

Let us finally mention that the coefficient b has a discontinuity for the aspect ratio £/L = 1: in this 
limit the corners disappear in the periodic geometry, and become four tt/2 corners in the open geometry. 
We refer to [29, 30] for a study of the corresponding universal terms. In the following, we shall focus solely 
on the case < l/L < 1. 
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2.2. Logarithmic terms 



CFT derivation. — The universal scaling function filjV) can be obtained as follows. We use complex 
coordinates w = x + ir on the cylinder (x £ [0, L)). Let C be a contour around the slit, and Wq = x Q + ir a 
point on the slit. Wc consider the infinitesimal transformation w > w + 5£Q(x — xq) inside the contour, and 
w i y zu outside. Such a transformation stretches the slit by 81 in the horizontal direction, while keeping L 
unchanged. The variation of the free energy is encoded in the T xx component of the stress-energy tensor as 

81 ( 

(T X x) dw (16) 



SF = 



2tt 



In complex coordinates we have T xx = T(w) + T(w). The integral can be obtained by mapping back the 
slittcd cylinder to the upper half plane H = {z £ C,Imz > 0}, where (T(z))g = due to translational 
invariance. The transformation law of the stress tensor is 



T(z) = K(z)) 2 TH + ^{ W (z),z}, 
where {w(z), z} denotes the Schwarzian derivative of w(z) 



{w(z),z} 



w"'{z) 



w"{z) 



(17) 



(18) 



w'(z) 2 \w'(z) 

The contributions of the holomorphic and anti-holomorphic parts of the stress-tensor are identical, and we 
get 

SF _ c 
~Sl ~ ~Vht 



S(z)dz , S(z) = (w'(z)) 1 {w(z),z} 
The conformal transformation from the upper-half plane is given by§ 



(19) 



w(z) 



2Z7T 



log 



1 + z z e 



2 p iirt/L 



z^- + e' 

and the integral (19) can be evaluated using the Residue theorem. We find 

SF _ TTC 

11 ~ 8Ltan(7r£/L)' 
which after integration yields the term of order one in (14), 

jWZ) = |log (sin^ 

Adding the additional contribution coming from the Cardy-Peschel formula (15), we finally get 



£ p - ail = - log 



L . 




— sin 
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+ 0(1), 



(20) 



(21) 



(22) 



(23) 



(v) 

where ajf is a non universal constant. Notice that such a "chord-length" scaling is quite common, as it 
appears for example in the study of the entanglement entropy [31], albeit with a different prefactor. Also in 
our case this is a subleading term, due to the presence of a line free energy. 

The open geometry can be treated as well. The conformal transformation from the upper half-plane is 
modified toll 



. . 2L 
w(z) — — arctan 

7T 



tan 



tt£ 
2L 



(24) 



§ A possible way to find this transformation is to map the cylinder to the plane, using u = e 2l7rw / L . The image of the slit is 
a circular arc with radius 1. This arc can be mapped to R + using a Moebius transformation f = (u — u±) / (u — U2), where u\ 
and U2 are the images of the two endpoints of the slit. The resulting geometry is sent to the upper-half plane using z = y/Q. 
This finally gives us z(w), which may be inverted to get w(z). 

|| The strip is mapped to the upper half-plane using u = e lnw / L , The image of the slit is again a circular arc with radius 1, 
which is mapped to a slit [0; i] using a Moebius transformation C(w), and mapped to the upper half plane using z = \J\ + £ 2 . 
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and we find 

+ 4 0) + o(l) (25) 

after a similar calculation. It is important that the £ up spins stand at the beginning of the open chain. 
If this is not the case the slit does not touch the left boundary any more, and the corresponding space- 
time geometry has the topology of an annulus, as opposed to the upper-half plane in our derivation. Global 
conformal invariance would then not be sufficiently strong to constrain the finite-size scaling function f(£/L), 
and this would make its determination much more difficult. 



£ p - a x l = -— log 



4Ltan»(^; 



Numerical checks. — We provide here some numerical evidence to support our main results (23) and (25), 
focusing on the Ising chain in transverse field. As is well known, such a chain can be studied by a mapping 
onto free fermions [32]. The calculation of the EFP essentially boils down to a determinant, as follows from 
the Wick theorem. See e.g Refs. [23, 29, 33] for the details. In particular we have 



det 

l<ij<i 



in the periodic case, and 



V = det 

i<i,j<e 



1 j 
2 



2 
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2L 



n(i-j + 1/2) 



1 



4L + 2 



CSC 



7T(t ~ j + 1/2) 

2L + 1 



4L 



CSC 



tt(x + j - 1/2) 
2L + 1 



in the open case. Here esc denotes the cosecant function 

1 

CSC X = — — . 



SIM 



(26) 



(27) 



(28) 



It is then straightforward to compute these determinant numerically, using standard Linear Algebra routines. 
The results are shown in Fig. 3, where we plot £ — a\l for system sizes up to L = 4096 and several aspect 
ratios t/L. As can be seen from the plot, the numerics agree very well with the CFT result, and even 




(a) Periodic geometry 



» -0.3 




(b) Open geometry 



Figure 3: Numerical extraction of the logarithmic terms (Eqs. (23, 25)). We show the data for system 
sizes L = 64,256,1024,4096 and the aspect ratios t/L = 1/16, 2/16, 3/16, 15/16. The regular spacings 
between the curves reflect the global scaling of the EFP, (c/8) logL for (a) and — (c/16) logL for (b). 



more so for the largest system sizes. Notice that in the periodic geometry, the invariance under I — > L — t 
is only exact in the asymptotic L — > oo limit. Indeed, we shall see in the next subsection that the first 
subleading correction is actually antisymmetric under £ — > L — l. Let us finally mention that we also checked 
with success our CFT formulas away from the ICTF (h = 1, v ^ 1). This provides strong evidence for the 
universality of the EFP. 
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2.3. An unusual subleading correction 

Let us now discuss lattice effects, and the possible subleading corrections they generate. In general these 
may be understood by adding perturbations by local operators to the CFT action. The operators can live 
in the bulk, or be localized near the boundary. In any case they need to be consistent with the symmetries 
of the physical system, and also be irrelevant (or at worst, marginal). Indeed, any relevant operator would 
drive the system to a different bulk CFT, or a different conformal boundary condition. 

Here we focus on a generic perturbation by the stress-tensor itself along the boundary, following [25] . 
The stress-tensor is arguably always one of the leading irrelevant operators, and we shall see that it actually 
produces the leading finite-size effects in our problem. Such a perturbation also admits an interesting physical 
interpretation [34]. Let us modify the CFT action to 

S^S+^ f dxT xx , (29) 

where the perturbation is localized near the slit. The coupling £ has the dimension of a length, and is 
typically of the order of one lattice spacing. Such a perturbation effectively pushes the conformal boundary 
to a distance £, to make up for the fact that the physical (lattice) boundary and the (RG invariant) conformal 
boundary have no reason to lie at the same exact position. For this reason, £ is nothing but the "extrapolation 
length" , a standard concept in boundary critical phenomena [24] . Let us also mention that this perturbation 
plays a key role in the study of quantum quenches [35] as well as the entanglement spectrum in fractional 
quantum Hall states [34]. 



CFT derivation. — Here we reproduce the main arguments of Ref . [25] , and refer to it for the detail. In the 
presence of a corner with angle 27r, perturbing by the stress tensor produces to leading order an interesting 
correction of the form log€, whose prefactor may be exactly determined for simple enough geometries. 
In general the contribution of the perturbation (29) to the free energy is given by the following power series 

Ai ? = -V'(— ^) — ! dxi [ dx 2 ... / dx n (T xx (xi)T xx (x 2 ) ■ ■ .T xx (x n )}. (30) 

„_i \ Z7r / n - Jslit ./slit J slit 



»W = i47T- ( 31 ) 



Let us explain the physical origin of the logarithmic correction on the simplest example of the infinite 
geometry. In this limit [L — > oo), the conformal transformation simplifies into 

£ z 2 -1 
2 z 2 + 1 ' 

The slit is the interval [— £/2,£/2] and its image through the inverse conformal mapping is M.+. Using 
T xx = T(w) + T(w): 

AF = -^ (T(w))dw=4r S(z)dz. (32) 



12tt 



w- 1 (-l/2+e) 



We have used once again the mapping to the upper-half plane and (T(z)) = 0. The regulator e in Eq. (32) 
is of the order of a lattice spacing, and acts as an ultraviolet cutoff. S(z) is here given by 

S{*) = l(^-Z + ^] (33) 



21 \2z 3 z 2 r 

For £ large we also have w~ 1 (—£/2 + e) ~ and w~ 1 (£/2 ~ e) ~ \f£jz, anc ^ tms auows to get the 

subleading contributions to the free energy. The term proportional to z in S(z) generates a correction ~ 1/e 
of order one, and the term proportional to 1/z 3 generates a correction in ~ e£~ 2 . The most interesting 
contribution comes from the 1/z term 

which is the £ 1 \og£ advertised earlier. Notice that (34) also yields a contribution proportional to £ 1 loge. 
This term is much less interesting to us, but it is important to take it into account for numerical purposes. 
The f _1 log^ is specific to corners with angle 2tt. It can be shown that for an angle 6, the corresponding 



9 



contribution to S(z) takes the form z 1 e / 7T , which can only produce a logarithm after integration when 
6 = 2tt. 

In the more general finite-size case we use the result of [25] . The contribution of interest is given by 
24tt \ ^ 



r Res [S{z);z 



x log I 



S{z) = («/(«)) 1 {w(z),z},(35) 



\w"(z c )\ 

where the sum runs over all the 2ir corners, their positions being w(z c ). w{z) is the conformal transformation 
from the upper half plane. The open geometry has only one such corner, and applying this formula we obtain 

cos(^) lo K £ 



9o{i/L) 



32 sin 



A similar calculation, with two corners, yields 
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x cot 



log^ 



(36) 



(37) 



in the periodic geometry. This term reproduces Eq. (34) in the limit t/L — > 0. Eqs. (37, 36) can be seen as 
"semi-universal" . While they are proportional to the extrapolation length which is a non universal quantity, 
the central charge also appears and the shape as a function of t/L is purely geometric: it does not depend 
on the specifics of the model. Another way to see this is to replace the dominant length scales L and I by 
L + e and I + e, where e and e are short distance cutoffs, of the order of a lattice spacing. Performing this 
substitution in Eq. (14) and expanding again, it is easy to see that the L~ x log£ terms remain unaffected. 



Numerical checks. — The extrapolation length is a so far unknown function of v on the Ising critical line. 
It should a priori depend on the boundary condition, but it has been found that £ = 1/2 for free boundary 
conditions in the ICTF [25]. To make contact with our predictions we consider the combination 



£ p — a\l - 



log 







— sin 




7T 





(38) 



in the periodic geometry, and a similar combination in the open case. We then fit this to do + d\L~ x log£ + 
diL~ x for various system sizes up to L — 4096 and several aspect ratios t/L. The results are summarized in 
Fig. 4, and show again a very good agreement with our prediction. We emphasize the importance of adding 
a last term diL~ x to remove the otherwise large finite-size effects. Finally, we also checked the validity of 
our result away from the ICTF, where the guess £ = \/{2v) is compatible with the numerics. An analytical 
argument for this value will be given in Sec. 3. 





ICTF. periodic' • 
CFT, periodic 































ICTF. 


open 




CFT, 


open 






t/L 

(a) Periodic geometry 



t/L 

(b) Open geometry 



Figure 4: Numerical extraction of the L~ 1 log£ term for the ICTF in both geometries, for aspect ratios 
t/L = 1/16, 2/16, . . . , 15/16. The agreement with the CFT predictions of Eqs. (37, 36) is close to perfect, 
with a relative error smaller than to one percent. 
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3. Relation to the theory of Toeplitz determinants 



E v = aj+- log£ + - |- -f- + 0{l/l), (39) 
£ =a x l-~ log* + a ( o) + + 0(1/*). (40) 



We now turn our attention to the limit L/t — > oo, so that £ becomes the only remaining length scale in the 
problem. The results of the previous section simplify into 

c . , , (p) £c log* 

- log £ + aX 

C , (o) . (,C log* 

— log t + a n H 

16 8 32tt * 

For the Ising critical line c = 1/2. £ is a so far analytically unknown function of ;/ on this line, but 
£(v = 1) = 1/2 for the ICTF. For simplicity we mostly focus on the ICTF point, but the generalization to 
the Ising critical line is straightforward. Interestingly, the determinants of Eqs. (26) and (27) now take an 
even simpler form: 

?v= n< Aet Oh-,-), (41) 

0<i,j<£— 1 

Vo= n *&£, A9i-j + 9i+j+t) , (42) 

0<2,j<e— 1 

with 

■ 9fc -^ + 27 r(fc+l/2)- (43) 

The factor (— l) k doesn't change the determinant, but has been added to make contact with the existing 
literature. An important feature of these determinants is that the matrix elements do not depend^ on 
the size *, and this allows for a more rigorous treatment, using known mathematical results. In the 
following we summarize some of these, regarding the asymptotics of Toeplitz (41) and Toeplitz+Hankel 
(42) determinants. This presentation is heuristic, and does not make any attempt at mathematical rigor. 
We refer to Refs. [36, 37, 38] for reviews. We then use these results to find the exact asymptotic expansion 
of the determinants corresponding to Eq. (43). Let us finally mention that the first three terms (but not the 
fourth) in Eq. (39) have already been derived by Franchini and Abanov [23] ; we will come back to this point 
later. 

3.1. The strong Szego limit theorem and the Fisher- Hartwig conjecture 

It is customary to see the matrix elements as Fourier coefficients of a certain generating function, usually 
called symbol: 

g k = [g]k = ^ r , 5(0) = E ^ ( 44 ) 

The following central result is due to Szego [39]. If g(<j>) is a sufficiently smooth, non vanishing function with 
winding number 0, the asymptotic expansion of detTi(g) — deto<i.j<£-i ([g]i-j) when £ — > oo is given by 

det(TH 5 )) ~ (G[g]f E[g], (45) 
where G[g] is the geometric mean of g(4>) 

G[g] =exp([log 5 ] ), (46) 
and E[g] is given by 

E[g] = exp ^ fc[log 5 ] fe [log 5 ]_ fe j . (47) 

The leading exponential term can be intuitively understood by noticing that Tt is almost invariant with 
respect to translations. Writing down Tg in Fourier space yields for large * 

[/|^C/,«diag[ 5 (0), 5 (27r/*),..., 9 (27r(*-l)/*)] , (Ufa = e" i2 ^ ; ( 48 ) 

1 Contrary to the determinants (26, 27), because of the fixed non-zero tjL aspect ratio. 
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provided the [g]k decay sufficiently fast. The determinant is then dominated by the geometric mean to 
the £— th power. Interestingly this result is exponentially accurate. There are however many interesting 
physical problems for which the [g]k decay slowly, and the symbol has singularities. This especially tends 
to happen when studying critical systems. Among them one can mention the famous Ising magnetization 
correlator at the critical point [21], the monomer correlator for dimers on the square lattice [40], full counting 
statistics [41, 42], entanglement in critical Id spin chains [45, 46, 47], just to name a few. The following 
result [48, 49, 50] will be useful for us. Consider generating functions that only possess pointwise "phase" 
discontinuities, parametrized as 

R 

9{<t>) = f{4>) I| exp (ip r arg 
r=l 

with — 7r < fa < fa < . . . < <f>R < 7r, |ft| < 1/2 and f((f>) sufficiently smooth, e.g. satisfying the hypothesis 
of the strong Szego limit theorem. We use the convention arg 2 6 (— tt;tt]. It is also possible to account for 
root-type singularities of the form (2 — 2 cos </>)", but we discard them here to make the discussion simpler. 
The asymptotic formula for deto<ij<^_i {[g]i-j) is modified to 

R 

det(T e (g)) ~ {G[f]ffE\g\ , c = -]Tft 2 (50) 

r=l 

where E[g] is a somewhat complicated constant. The most interesting physical piece in such formulae are 
the power law correction £ u , because they usually correspond to universal physical terms. 

Such a prescription is known as the Fisher-Hartwig conjecture [22], and many instances of this 
conjecture and various generalizations are by now theorems [51, 52]. In the following we will also 
need similar results for the closely related Toeplitz+Hankel determinants of the form det((T + H)i) — 



(49) 



3.2. Rigorous asymptotic expansion of the EFP 
For the Ising critical line the symbol reads [23] 



r±s 1 1 / • 

1\(f>) = - + — exp [ —1 arctan 



v tan — 
2 



(51) 



and at the ICTF point (y = 1) it simplifies into 

fM) = \ + SSn(CO f /2)) e-^. (52) 

It has an obvious phase discontinuity, due to the change of sign in cos(</>/2). It can be parametrized as in 
Eq. (49), with R = 1, fa = and ft = P = -1/4. f(fa is given by 

fW = IsWI = ^2 + 2|cos(0/2)| (53) 

and has winding number zero, as it should. Armed with this, it is now straightforward to apply the results 
of [50] , and get the leading terms in the asymptotic expansion of the emptiness formation probability. 

Periodic case. — We recover here the result of [23] , which reads 

£ = ai £+^log£ + a^ +o(l), (54) 
compatible with the CFT prediction. The non universal constant is given by 



a o P) = -5>([log/] fe ) 2 -l°g(G(3/4)G(5/4)) = 0.0951599254(1) (55) 

fe=i 



G is the Barnes G-function [53], solution of the functional relation G(z + 1) = T(z)G(z), where T is the 
Euler Gamma function. The numerical value has been obtained by brute force numerical evaluation of the 
first 5000 terms in the sum. The logarithmic term agrees with CFT, and it is easy to check that it remains 
unaffected by a change in v. 
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Open case. — We apply here the result of Ref. [50], and refer to it for the general formula. Let us just 

j_ 

32 

then find 

f. = a-if- 

32 

with the non universal constant given by 



mention that the critical power law exponent is modified to u> — — ^ — |/3 2 = ^ in this particular case. We 
£ = ai^T^log^ + 4 o) +o(l), (56) 



1 oo „ 



0) 

o 2 ^ » 4n 

k=l 



G(3/4) 2 G(5/4) 

G(l/2)2 7 /32 7r l/8 



-0.02149124775(5). (57) 



C is the Catalan constant. Once again, this result agrees with CFT. As in the periodic case, the logarithmic 
term does not change when varying v. 

3.3. Possible subleading corrections 

Let us now discuss possible sources of subleading corrections to the general Fisher-Hartwig formula. If one 
relaxes the condition on the (3 r in Eq. (49), the parametrization is not necessarily unique anymore, because 
each j3 r can be shifted by arbitrary integers /3 r — ¥ f3 r + n r , the zero winding condition on / being ensured by 
n r = 0. The following result, often dubbed "generalized Fisher-Hartwig" has been conjectured by Basor 
and Tracy [51] and proved recently by Deift, Its and Krasovsky [52]: 

/ R 

det(T e (g)) ~ (G[f\f £ ^M,{nr}) E[g . j K}) = _ ^ r + n r ) 2 (58) 

{n r } r=l 

where the sum runs over all possible n r € Z subject to the condition J^^Li n r ~ 0. It is important to stress 
that the symbol ~ stands for equivalent. This formula therefore only differs from Eq. (50) when several 
different representations have the same exponent Rcw({(3 r }, {n r }). 

It is however tempting to assume that all the terms in (58) appear as subleading corrections to the 
determinant. Several studies [54, 55, 56, 41, 42] supported this claim, by showing that at least the first 
subleading branches can be identified in the asymptotic expansion. Moreover, each branch can separately 
be improved by adding corrections in the form of a power series in I . Based on explicit computations of 
the first subleading terms, Kozlowski conjectured [26] that this prescription yields, provided f(<p) is smooth, 
a full asymptotic expansion of the determinant. More precisely, it is given by 

det(TH 5 )) = (G[f}f £ £»(iM,{n r} ) E[g . {nr}] f 1 + g a f 0r} {nr} lA . (59) 

{n r } \ i=l / 

In our case the generating function has only one Fisher-Hartwig singularity, which means there is no other 
subleading representation+ . Because of the cusp at <j> — n, /(</>) is continuous but not smooth, hence our 
i^ 1 log£ does not contradict Eq. (59). The Riemann-Hilbert analysis of Ref. [26] can however be generalized 
to account for less regular symbols, by introducing the parametrization 

5 (0) = /(0)(l + z)-^ ) (l + l/z)^ ) , z = e», (60) 

where v and v are analytic functions of z in some neighborhood of the unit circle, and chosen in such a 
way that /(</>) be also analytic*. For symbols g{4>) that can be written in such a way, the Riemann-Hilbert 

(i) 

approach allows to obtain the form of the subleading corrections: in general each coefficient a \p } { n } m 
Eq. (59) should likely be replaced by a polynomial of degree at most i in log£ One can check that writing 



1/(2)= -p-J^Tbiz + lY (61) 



p=l 



+ We point out a slight inaccuracy in Ref. [23], where the authors considered a parametrization with /3i + 1 and conjectured 
a subleading term in £-(/3i+i) 2 +« = £-1/2 in 

the LEFP. Such a parametrization is however forbidden, because /(</>) doesn't 

have winding number zero in this case. 

* I am grateful to Karol Kozlowski for explaining this to me. 
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OO 
33=1 



(62) 



with f3 — — 1/4 as before preserves the phase discontinuity at 9 = tt. Notice that ry p = for all p would give 
back the previous parametrization of Eq. (49). The rj p have to be set so as to smoothen the more regular 
part of the symbol f(<j>) around the singularity at 6 — ir. Let us determine the first few for the Ising critical 
line of the XY chain, where the symbol is given by Eq. (51). 771 = — j- is obtained by imposing f'(n) = 0, 
and f]2 = Tli/ 2 is the only value that makes the second derivative f"{ii) finite. It turns out the knowledge of 
■qi is sufficient for our purpose. Pushing further the analysis of Ref. [26] we get a term of the form 

1ok£ 



(63) 



(64) 



This yields an extrapolation length 
1 

consistent with our previous determination at v = 1. We also study a generalization, relevant to the full 
counting statistics, in Sec. 5.2. Finally, let us mention that a similar structure of logarithmic corrections 
has already been studied for closely related Fredholm determinants of the generalized sine kernel (see e.g. 
Refs. [43, 44]). 

Since our results are based on heuristic arguments, we are unable to prove the existence of the £~ x \og£ 
term. To numerically resolve this issue, it is convenient to form the combination 

E(£) = (S - - b log£ - a ) X £ (65) 

either in the periodic or open case. Provided the CFT result is correct, this quantity should scale as 

E{1) =7logi + 5 + o(l). (66) 

Fig. 5 shows some numerical data for the ICTF in the periodic case, where E(£) is plotted against log£, 
for various system sizes between £ = 128 and £ = 16384. The data follows an almost straight line, and 
confirms the scaling of Eq. (66). Our best numerical estimate is 7 = —0.00994(2) which agrees with the 
predicted 7 = — l/(327r) within less than 0.1 percents of relative error. We also observed that this estimate 



-0.06 - 



-0.07 



-0.08 



-0.09 



-0.11 - 



Numerics )K 
Best fit: -0.009934 x x - 0.009702 




log<? 

Figure 5: E(£) for system sizes £ = 128,256,512,1024,2048,4096,8192,16384. 
function of log£ and agrees very well with the scaling predicted from CFT. 



The data is plotted as a 



can be improved by adding further subleading corrections. We found that adding two terms of the form 
b_ 2 £ 2 log 2 £ and a_2^~ 2 reproduces the numerical data very well, and brings the relative error on 7 under 
4.10~ 6 . A similar (but not identical) scaling has also been observed [57, 25] for the fidelity considered in [58], 



14 




where a closed-form formula can be derived and a systematic asymptotic expansion be performed. In this 
case there is a l~ x \og£ as well, but the subleading corrections take the somewhat simpler form of £~ v \og£ 
and £~ p terms, with p > 1. 

Let us now summarize our findings. In general we expect the LEFP to scale as 

a P ,k (log £) k \ XT? (67) 

fc=0 / 

where the prefactor of the t~ v is a polynomial of degree at most p in log t. Such a scaling is natural in terms of 
higher-order corrections to the Toeplitz determinant, but also in CFT. Indeed, taking into account the other 
contributions of the perturbed expansion (30) leads to a similar structure for the corrections, although the 
precise determination of all terms becomes very quickly cumbersome. Note however that not all logarithmic 
terms in the expansion (67) are semi-universal, in the sense defined in Sec. 2.3: making the substitution 
£ — > £ + e in Eq. (67) and expanding again, we observe that a (log^) fc £~ p remains unchanged only if all the 
lesser subleading corrections are of the form (\og£) k £~ p , with k' < k and p' < p. In particular this implies 
that the (\og£) k £~ k are automatically semi-universal, but the other terms usually not; the former would be 
the most interesting to study within CFT. 

4. The ferroelectric string as a boundary that breaks scale-invariance 

We now turn our attention to the XXZ chain, and related models. The Hamiltonian reads 

H =Y, (°Ku + a j a j+i + A ^J + i) • (68) 

3 

It is well established that the low energy physics of such systems is described by a Luttinger liquid CFT in 
the range — 1 < A < 1. In imaginary time this is a Gaussian model, with action 

S = ^- I (V/i) 2 dxdy , h ee h + 2irr (69) 



47T w 

The height field h lives on a circle of radius r. The physical parameter that governs the decay of the 
correlations is the compactification radius R = ry2K. For the XXZ chain it is given by 



R{A) = ^ 2 - ~ arccos A. (70) 

Our analysis will focus for technical reasons on the free fermion point A = 0, although we expect our 
reasoning to apply to the whole critical region. We also assume half-filling throughout the entire section. 

4--1- Scaling of the EFP 

Many exact results are available for the scaling of the EFP, both for the infinite and semi-infinite geometries. 
Generically we have 

£ = a 2 £ 2 + ai £ + b a log£ + o(l), (71) 

so that the leading term is proportional to £ 2 , and the analysis of Sec. 2 does not apply anymore. We 
explain the precise physical origin of this leading term in the next subsection. Intuitively the difference with 
Ising comes from the fact that the ferroelectric string injects a lot of particles (or up spins) in a system 
where the total number of particles is fixed (or that has zero total magnetization) . We also want to provide 
numerical evidence for the universality of the subleading logarithms in 4.3. To try and achieve both goals it 
is convenient to introduce another free fermionic spin chain 

H d = ^ [°h ( C2 J- 2 + + 3C 2? + ° 2 3 + 1 + C 2j + 2) + 4.7+1 ( C2 J + C 2j + 1 + c 2j+2) +fl.c(72) 

3 

Just like the XXZ chain has the same eigenvectors as the transfer matrix of the six-vertex model, H has the 
same eigenvectors as the transfer matrix of the dimer model on the square lattice [59, 60, 61]. Hence we dub 
(72) the "dimer chain" . In terms of dimers a f spin (or a fermion) is either an even vertical link occupied by 



15 



a dimer, or an empty odd vertical link. The ferroelectric string in the spin language corresponds here to an 
alternating sequence 1 1010 . . . 10) of links occupied (1) by a dimer, and empty links (0). The EFP of both 
periodic chains can easily be calculated as a determinant V = deto<ij<i-i(mij) using Wick's theorem [14]. 
The matrix elements are given by 



sm 



\ii-3) 



L sin J (i 



.i) 



(73) 



for the periodic XX chain at half-filling. In the limit L oo the corresponding symbol takes the form 



i 



sgn (cos < 



7T < 4> < 3tt/2 
otherwise 



(74) 



so that it vanishes on two whole intervals of [0;27r]. This is a big difference with the Ising case, as we 
cannot apply Fisher-Hartwig related results to this determinant. Here the enhanced scaling proportional to 
I 2 follows from Ref. [62]. In the dimer chain the matrix elements are given by 

fa cob[§(1 + (-1V-'T 



2 



L 



9SQ 



cos z q 



■coB[q(i-j)] 



(75) 



Here Q = { —tt + (2m + 1)tt/L \ m = 0, . . . , L/2 }. This matrix is only block Toeplitz, because of the factor 
(— l) lJ . There arc also similar formulae for the open chains. Finally, note that such determinants are closely 
related to the so-called "gap-probability" encountered in the study of random matrices [63] . 

4-. 2. Imaginary time picture, and the arctic circle phenomenon 

Before presenting our results, let us mention that it is possible to look at the imaginary time picture for 
a weakly magnetized string, where bosonization applies [17]. It was shown that already in this limit, the 
slit shown in Sec. 2 becomes an ellipse, where the degrees of freedom are frozen. The LEFP then becomes 
proportional to £ 2 . It was however noted that this picture agrees only qualitatively with the numerical data 
at maximum magnetization. Here we determine the region numerically for finite but large I. 

We first recall the height mapping for dimers. Let us associate an integer height with each plaquette of 
the lattice, as illustrated in Fig. 6. The heights are defined as follows. Set the value on the lower leftmost 
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(a) 6 x 6 square lattice. 



(b) Emptiness formation probability 



Figure 6: Height mapping, (a): 6x6 square lattice with open boundaries, (b): a sequence |tttt) = 1 1010) 
is imposed at the left of the fourth row in a semi-infinite grid. The dimers are the red thick lines, the empty 
links are the red zigzag lines. Due to the hardcore constraint, the occupation probabilities of several other 
links are set to p = or p — 1 (in blue). The frozen region corresponds to a maximal slope for the height. 



face to be zero. Then, turning counterclockwise around a site of the even (resp. odd) sublattice represented 
by black squares (resp. white squares), the height picks +3 (resp. —3) when crossing a dimer, —1 (resp. +1) 
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otherwise. It is easy to see that the flat configurations should dominate the others: after coarse-graining 
the long distance physics is well described by the gaussian action of Eq. (69). In contrast, the ferromagnetic 
string (shown in red in Fig. 6b) which we impose corresponds to a maximum slope on the height field. This 
has far reaching consequences, due to the hardcore constraint: Fig. 6b shows (in blue) all the links whose 
occupation probabilities are automatically set to and 1 for a finite slit length I = 4. The dimers are 
therefore completely frozen, and the height field has a maximal slope in the whole frozen region. 

In the limit £ — > oo such constraints become highly non trivial, as the frozen region actually extends 
further than one could naively think. Fig. 7 shows the frozen region in the infinite and semi-infinite 
geometries, for a slit length £ = 56. The pictures are obtained as follows. First we calculate the probabilities 
associated to dimcr occupancies on all links. The procedure to get them is described in Appendix A, and 
makes extensive use of the Pfaffian solution of the dimer model. Second we associate a mean height to each 
plaquette using these probabilities)!. Third we define a discrete squared gradient 

(V d /i) 2 = [h(i x ,i y + 1) - h(i x ,i y - l)] 2 + [h(i x + l,i y ) - h(i x - l,i v )] . (76) 

The maximum gradient is (Vd/i) 2 = 16. We draw in blue the frozen region where the gradient is maximal 
and in black the regions with zero gradient. Intermediate regions are shown as a mixture of the two colors. 
Fig. 7a (resp. 7b) shows the results for the frozen region in the infinite (resp. semi-infinite) geometry. We 
also show in Fig. 7c an example of height profile. 




(a) frozen region (b) frozen region (c) height profile 



Figure 7: Arctic phenomenon generated by the ferroelectric configuration ft • ■ ■ t- ( a ) : Slit of length £ = 56 
in the infinite geometry. The data is shown on a 112 x 84 rectangle, and the frozen region is in blue, (b): 
Slit of length £ — 56 in the semi-infinite geometry. The data is shown on a 112 x 168 rectangle, (c): Slit of 
length £ = 56 in the semi- infinite geometry. Same setup as for (b), with the high heights shown in red. 

The limiting shape is a non trivial curve with several cusps. Notice that outside of the frozen region 
the dimer occupancies are still constrained (i.e p ^ 1/4), although they go to this value sufficiently far 
away. Using this we can interpret the leading term in the EFP: it is simply due to a deficit in free energy 
coming from the arctic region, as well its proximity, where the free energy is not constant with the position. 
Interestingly, the limiting shape for the semi-infinite geometry with slit length £ seems to match perfectly 
half that of the infinite geometry with slit 2£. This observation is consistent with the exact results of [19, 15], 
which give 0,2 (semi — inf) = 2a2(inf) in the XX chain. 

Inspired by the results of Sec. 2, it is tempting to look at the arctic regions for a finite aspect ratio 
£/L, to see if the frozen region is modified. We present in Fig. 8 some numerical results regarding this 
question in the open geometry (the periodic geometry is similar). We observe that, as long as £/L remains 
relatively small (figs. 8a and 8b), the picture looks more or less the same, up to a global dilatation. However 
the leading coefficient seems to change slightly (anticipating on table. 2), when varying the aspect ratio, so 
that this observation is possibly only approximate. Also the limiting case £/L — 1/2 is interesting, because 

(J For example the height picks ±(3xp-lx(l- p)) = ±(4p — 1) when crossing a link occupied by a dimer with probability 
p, and empty with probability I — p. 
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(a) I = 16 



(b) I = 24 



(c) / = 32 



Figure 8: Arctic phenomenon for an open system of total length L = 64. The frozen region is shown for 
three different slit lengths. 

the sequence of L/2 consecutive up spins has to be followed by L/2 consecutive down spins, due to the 
zero magnetization constraint. In this case the frozen region even cuts the imaginary time picture in two 
disconnected parts (see figure. 8c). 

Such a behavior is yet another manifestation of the arctic phenomenon, first discovered in the study of 
dimers on the Aztec diamond [64]. Compared to the standard rectangular geometry, the dimers are highly 
constrained near the corners of the diamond, and these constraints propagate on macroscopic scales. In 
terms of heights the slope is maximal on all sides of the diamond, and the authors of [64] managed to prove 
that all degrees of freedom outside of the arctic circle were frozen. An illustration of this phenomenon is 
shown in figure. 9a. Note that the Aztec diamond can easily be generated by imposing maximally magnetized 




(a) Arctic circle on the Aztec diamond (b) Aztec diamond and emptiness formation probability 

Figure 9: Arctic circle phenomenon, (a): Frozen region on a 32 x 32 Aztec diamond. All dimer occupancies 
are set with probability one outside a circle in the thermodynamic limit, (b): The Aztec diamond (thick 
black lines) can also be obtained by imposing the boundary conditions 1rtt1rtt44i-UU44- & t the left and right 
of a I x I square (in red). In terms of dimer occupancies the boundary condition reads 1010101001010101. 
Inside the circle lives a field-theory, whose exact properties are so far unknown. 

boundary conditions on a square: in terms of transfer matrix the partition function for dimers on the Aztec 
diamond of diameter L x reads 



provided L y > L x . T is the transfer matrix of the dimer model, acting on the vector space of dimension 
2 L y generated by the dimer occupancies along a vertical line. See figure. 9b for an illustration in the case 



z = (t... U •••4 |T £ -|t...t4 ■■■i), 



(77) 
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L x = L y = 16. Let us finally mention that the arctic curve separating the frozen region from the fluctuating 
(CFT) region can usually be understood in terms of a entropy minimization principle [65, 66]. It would be 
interesting to apply this method to our problem, and attempt an exact calculation of the regions shown in 
figures 7a, 7b. 



4-- 3. Logarithmic terms 

As for the conformal boundary case, the logarithmic terms in the EFP are potentially the most interesting, 
as it is reasonable to assume that they might be universal. Many exact results are already available, they 
are summarized in table. 1 for the XXZ chain. We also show the terms for the six-vertex analog of dimcrs 
on the Aztec diamond, namely the six- vertex model with domain wall boundary condition [67, 68, 69]. The 



A 


6(A), periodic geometry 


6(A), open geometry 


Arctic "circle" 





1/4 


1/8 





1/2 


5/36 


-1/72 


-1/18 


general 


i (2R 2 + &/R 2 - 7) 


? 





Table 1: Summary of the logarithmic terms in the XXZ chain. The general result for 6(A) in the periodic 
geometry has been conjectured and numerically checked in [17]. The others are exact [14, 15, 19, 70, 71]. 



prefactor of the logarithm is always given by a simple function of the compactification radius, and does not 
depend on non universal features such as the Fermi-speed, etc. This is an additional hint for universality, but 
unfortunately the field-theoretical analysis of the exterior of the frozen region is not so easy to perform. The 
non frozen degrees of freedom should still have long-range correlations. The field theory describing them, 
although not strictly speaking conformal, should therefore still be massless. We speculate that the stiffness 
k of Eq. (69) might depend on the position, at the exterior of the frozen region. A possible way to show 
this would be to exploit some exact results for the dimers Green's function on the Aztec diamond [72], and 
determine this way the asymptotic behavior of correlation functions, in the spirit of Fisher and Stephenson 
[40]. This however falls outside the scope of the present paper. 

Such logarithmic terms can also be compared with those appearing on a rectangle, which are simpler 
from the CFT point of view. There is however an important subtlety that needs to be taken into account. 
Indeed using the height mapping, it is easy to see that the boundary heights are not exactly the same on the 
two sides touching a given corner. See figure. 6a, that shows a coarse grained height (h) = ±1/2 depending 
on the side. This small mismatch has a tremendous impact in the thermodynamic limit, and needs to be 
dealt with by adding a harmonic function in the action, that encodes the height shift [88]. In the language 
of boundary CFT, this corresponds to the insertion of 4 (magnetic vertex) boundary changing operators 
[89] . It is then straightforward to apply the general analysis of Ref . [90] , and one even gets this way the full 
partition function 

Zrect-W e/2 » T = l — , (78) 

[v(t)] Lx 
with central charge c = 1. Here v depends on the parity of L y : v = 3 for even L y and v = 4 for odd L y \]. 
At the free fermion point R — 1, and Eq. (78) reproduces the exact result of Ferdinand [92]. 6% and 9\ are 
the third and fourth Jacobi theta functions, given by 

3 (r) = 5>^, (79) 

feez 

ft The difference when L y is odd is related to the order in which the height shifts are applied: we have the boundary correlator 
{V ct (z 1 )V- a (z2)V c ,(z 3 )V- a (z4)) (resp. {Va(z 1 )V ct (z2)V- a (z 3 )V- a (z4)}) for L y even (resp. odd). Here z 1 (resp. z 2 , z 3 , 24) 
denotes the position of the bottom left (resp. top left, top right, bottom right) corner, and a 2 = R 2 /16 [88]. See [91] for a 
similar even-odd effect on the cylinder. 
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Wt) = 5](-l) fc e I7rfc2r , (80) 
and rj is the Dedekind eta function 

oo 

t](t) = e- MT/12 Y[ (1 e 2 " feT ) . (81) 
fe=i 

Notice the similarity of the logarithmic term (i? 2 /4— 1/4) logL in Eq. (78), with those summarized in table. 1. 
However, the rectangle is somewhat less rich, as it doesn't have electric vertex contributions proportional to 
l/R 2 . 

Here we pursue a more modest goal, by numerically studying, at the free fermions points, both the XX 
chain and the dimer chain, while also varying the aspect ratio x = i/L. The results for a periodic chain 
are shown in table. 2. Of course the non universal terms are not the same for the two chains. However, the 





XX chain 


dimer chain 


coef 


x = 1/8 


x = 1/4 


x = 1/3 


x = 1/2 


x = 1/8 


x = 1/4 


x = 1/3 


x = 1/2 


a 2 


0.3498456 


0.3604522 


0.3729935 


0.4262783 


0.3984028 


0.4139014 


0.4321317 


0.5047001 


bo 


0.2500058 


0.2500012 


0.2499928 


0.1666604 


0.2500008 


0.2499955 


0.2499783 


0.1666236 


a 


0.3468585 


0.3297074 


0.3065721 


0.2556067 


0.3295512 


0.3230345 


0.3161439 


0.3289082 



Table 2: Scaling of the LEFP in the periodic XX and dimer chains, for various aspect ratios x = t/L. The 
data is extracted via a fit to the form c^ 2 + bo log I + ao for system sizes between L = 80 and L = 96. 



logarithmic ones stand compatible with table 1, as long as the aspect ratio is smaller than 1/2. Precisely 
at x = 1/2 the logarithmic terms have a sudden jump, and go to a value that is very likely to be 1/6. We 
also performed a similar analysis in open chains, where we found logarithmic terms with prefactor 1/8 for 
x < 1/2 and 1/12 for x — 1/2. This is another evidence in favor of universality, as these terms seem to not 
depend so much on all the details of the arctic region. Their analytical determination is left as an important 
open problem. 



5. Some related problems 

We discuss in this section some possible extensions of this work, as well as some related questions. 



5.1. The other Ising critical line in the XY chain 



This line with 7 = 1 and h — — 1 also belongs to the Ising universality class, and has also been studied by 
Franchini and Abanov [23]. They showed, in the infinite geometry, that the LEFP scales as 

1 



£ p = ail 



\og£ 



Ap) 



l -l/2 



Hi ' " VI 

The slightly different scaling is mathematically explained by the fact that the symbol 
1 1 



exp 



- / arctan 



i^tan 



(82) 



(83) 



has an additional singularity at (f> = 0, and subleading representations have to be taken into account (see 
Eq. (58)). Our analysis can in principle be extended to explain this scaling. First, notice that the EFP in 
this chain (h = — 1) can be recast as the probability of observing a sequence 1 44- • • • X) i n the h = 1 chain. 
Then, we have 

|U...4),= 2^a | -(H)x- I (84) 
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= E (-i)^kW (85) 
{ CT ?=-h<-} 

± [free}*, (86) 
where iVj, counts the number of down spins. Since the free boundary condition is unstable under the RG flow, 
this should actually renormalize to a (stable) fixed boundary condition, and we have to introduce boundary 
changing operators. In the infinite geometry they do not affect the Cardy-Peschel term because the external 
boundaries are at infinity, but they influence a lot the subleading corrections. For example we still expect a 
log€ term to appear in (82), but with a prefactor modified to 

b. 1 =-£(c-16h {+ ), (87) 

where hf + — 1/16 is the dimension of the operator changing the boundary condition from free to +. We 
refer to [25] for the method. The changes are even more dramatic in the semi-infinite geometry. Applying 
the Cardy-Peschel formula (15) we get 

So = M + (4/if+ - ^) log^ + 4 o) + o(l), (88) 

so that the prefactor of the logarithmic term is modified to 7/32. To our knowledge this exponent cannot be 
obtained using known theorems on Toeplitz+Hankcl determinants. We also expect the presence of a € _1 \og£ 
correction among the subleading terms. 

5.2. Full counting statistics 

Another closely related problem is that of the counting statistics of magnetization in spin chains [73]. What 
makes this type of problem interesting is that it is motivated by recent experiments in cold-atom systems. 
In particular, it is sometimes possible to use the quantum noise to measure the full probability distribution 
of fluctuating observables [74, 75]. 

For exact studies it is usually more convenient to study the magnetization in the basis of the eigenstates 
of the a z (see, however, Refs. [76, 77]). The statistics can be conveniently expressed in terms of the following 
generating function 

Jtt(A) = (e-^ES-fr-?)/*) = £ p m e-™\ (89) 

m>0 

where p m is the probability of having m spins down (out of £) in the z-basis. As we shall see, x^(A) is a 
natural generalization of the EFP. All cumulants are recovered as 



dno gxe (x) 



<^>e = W dXP 



(90) 

The generating function can be recast as a Toeplitz determinant x^(A) = det([g,\]i-j), with symbol [75, 78] 



, JN l + e~ A l-e~ A . . 
9\W = 2 1 2 6XP ' _zarctan 



j/tan — 
2 



(91) 



Next we apply the analysis of Sec. 3.3, limiting ourselves to e A 6 [— 1; +oo] for simplicity. The symbol can be 
put in the form (60), with f3 — — arctan[tanh(A/2)]/7r and r/i = — tanh X/(4irv). The asymptotic expansion 
reads 

- logx^(A) = ai{\)t + 6 (A) log* + a (A) + 6_i(A)r 1 log£ + ©(r 1 ), (92) 
and the coefficients are given by 

oi(A) =[log| 5A |]o (93) 

OO 

«o(A) = -EM[ lo slSA|] fc ) 2 -log(G(l + /?)G(l-/3)) (94) 
fc=i 



&o(A) = n arctan 2 



, A 

tanh — 
2 



(95) 



fc-i(A) = - tanh A (96) 

2ttv 
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ai(A) was studied over the unit circle |e A | = 1 in [78], and ao (A), bo(\) are extensions for v ^ 1 of the results 
of Ref. [75] for the ICTF [y = 1). The last term 6_i(A) is new. 

Two special values of the generating function are of interest, in addition to the obvious x^(0) = 1- First, 
lim^oo XiW is nothing but the emptiness formation probability, and we recover b = 1/16 = c/8, as well 
as 6_i = — l/(327w) = — £c/(87r). Second, one can show, using a Kramers-Wannier transformation, that [75] 

Xl&in) oc (afaf > ~ r 1 ' 4 (97) 

and one recovers the result of [82, 83]. In CFT language the exponent is twice the scaling dimension (1/8) of 
the Ising spin held. Note that it is possible to interpret the A dependence of &o(A) and 6-i(A) within CFT. 
In terms of the Majorana fermion fields, the generating function can be expressed as 

XtW = ^exp (^-i^J o dxil>(x)$(x^J \ . (98) 

where ip and if) are the analytic and anti-analytic components of the Majorana fermion. The argument of 
the exponential can be seen as a perturbation of the CFT action, and has been studied in Ref. [84]. The 
imaginary time pictures shown in Fig. 2 remain valid. There is however an important subtlety: the term 
Jq dx?p(x)'ip(x) is marginal along the slit, so that the conformal boundary condition are not free anymore. 
To determine them properly, one can fold the model about the slit [85] , and then bosonize the two resulting 
Ising copies. The corresponding boundary condition, for a given A > 0, is then part of the one-parameter 
family of continuous Dirichlet boundary condition investigated in Ref. [85]. We refer to to [84, 85] for a 
thorough discussion. 

From the generating function it is possible to reconstruct the cumulants. For example the hrst three in 
the ICTF are given by 

' ' \ (99) 



(to 2 ) 



2 7T 



h 1 

4 2tt 2 



log^ + cst + 0^" 1 ) (100) 



(to 3 ) c = ±1 + est - JL ^ + Or 1 ) (101) 

The linear scaling of all the cumulants means that the probability distribution p m is gaussian. Note that 6 
(resp. b-i) is symmetric (resp. antisymmetric) under to x — > x , hence log-^ (resp. i~ l \ogi) terms can 
only appear in even (resp. odd) order cumulants. 

It is useful to compare this to the XX chain case (see e.g [81, 41]). Here an important difference is that 
the magnetization is conserved, and the system is equivalent to free fermions without pairing. The symbol 
reads 

1 + e~ A 1 - e~ A 

9\W = 2 1 2 Sgn *- C ° S ^ ~~ C ° S kpS> ^ 102 " > 

and has two Fisher-Hartwig singularities at (f> = Uf and <f> = 2n — hp. kp is the Fermi momentum. The 
corresponding generating function is given by 

-log X£ (A) = — t+~logt + 0(l), (103) 

with subleading corrections taking the form (59). The result becomes singular when A — > oo; we have seen 
that in this limit — logx^ is proportional to £ 2 , not £. This is an important difference with the Ising case, as 
the counting statistics and the EFP have here a very different behavior. The cumulants are given by 

(m) c =— I (104) 

(m 2 ) c = -^log£ + (D(l) (105) 

The higher-order cumulants (m p ) c are either zero for p odd or 0(1) for p even. All these can also be 
understood from simpler bosonization arguments. Interestingly, the even order cumulants can also be used 
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to reproduce exactly the entanglement entropy [79, 80, 81], and therefore S = (tt 2 /3) (m 2 ) c at the leading 
order. Note that while this connection can be generalized to fermions with pairing (that describe the 
Ising chain or superconducting states), it is necessary to count the fermionic-quasiparticles appearing in the 
entanglement Hamiltonian [80, 81], not the original fermions as we did in (89). Hence the result we have 
established for the Ising chain does not seem to have any relation to the generating function whose even-order 
cumulants reproduce the entanglement entropy. This observation is consistent with the general results for 
the corrections to scaling [86, 87]. 



5.3. Classical Mutual information in a spin chain 

Let us consider the same bipartition as before for a spin-1/2 chain. Since Trp^ = 1, each diagonal element 
(erlp^lcr) of the reduced density matrix may be viewed as a probability. It is then natural to consider the 
classical Renyi-entropy of these probabilities 



Sn(L, I) = — - In V [(a | PA | a)] n , (106) 




from which the classical Shannon entropy is a limiting case n —> 1. Here the sum runs over all 2 possible 
spin configurations in subsystem A. This entropy is not to be confused with the Renyi entanglement entropy 
(REE), as the reduced density matrix is not diagonal. However when £ = L, it is nothing but the REE of 
two-dimensional Rokhsar-Kivelson states in a cylinder/strip geometry [61], and has been further studied for 
this reason [29, 30, 88]. For £ ^ L the correspondence doesn't apply anymore, but most of the techniques 
developed previously can still be applied. In particular, the general scenario of Refs. [29, 88], that predicts 
two different behaviors as a function of n separated by a phase transition, should still hold. This is interesting 
because in the Ising chain the limit n — ¥ oo of this entropy is nothing but the logarithmic emptiness formation 
probability. Before presenting some possible consequences of this connection, it is convenient to introduce 
the Renyi Mutual information 

I n = S n {L,£) + S n (L,L-£)-S n (L,L). (107) 

This combination filters out the line free energy contribution, so that the universal logarithms are now the 
leading terms. It was shown that for n > 1, all universal properties are dominated by the n — > oo limit [29]. 
Using Eq. (6) we predict 



L ir£ 



I n = - x log - sin — + est. + o(l) (108) 

8 n — 1 \7t L J 

for a periodic chain, and a similar result in the open chain. The limit n = 1 in the periodic case was studied 
numerically in Refs. [93, 94], and this is usually a challenging problem in CFT. However, we notice that the 
conformal scaling of Eq. (108) still seems to apply, albeit with a different and non trivial universal prefactor 
[93, 94] (see also Eq. (III. 181) in Ref. [57] for an even more precise determination in the limit L — > oo). 
Similar mysterious numbers have been observed numerically for the Shannon entropy of a full chain [29, 30]. 
It would also be interesting to study in details the subleading corrections at n > 1 and n = 1, to see if the 
unusual L^ 1 \ogL term is still present as a subleading correction. 

For the XXZ and dimer chains the limit n — > oo doesn't give the emptiness formation probability, but 
rather the probability of observing an antiferromagnetic string p(f4- ■ • • ti)- At the free fermions point the 
corresponding determinants are also simple and take the form p = deti<i.j<e (mij) with 

sin[f(i-j)] 



for the XX chain and 



mi 



St + H^> £g -H(. + (-i^)] 

2 L tt Q Vi + cos 2 ? 

for the dimer chain. Once again Q is the set Q = { —n + (2m + l)ir/L \ m = 0, . . . , L/2 }. Here we have only 
given the determinants corresponding to a periodic chain, but there are similar expressions in the open case. 
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It might be possible to study exactly the infinite limit L — > oo, but the situation is more complicated because 
the corresponding matrices are only block- Toeplitz. To the author's knowledge no mathematical result can 
be applied to these determinants as of now. It is however possible to perform a CFT analysis similar to the 
one presented in Sec. 2. For example the leading logarithmic terms can easily be predicted in both chains. 
Since the configuration with alternating spins is expected to renormalize to a Dirichlet boundary condition 
for the boson field, we predict a term in | log £ in the infinite geometry. In the semi-infinite geometry one has 

to take care of the height shift discussed in Sec. 4.3, and one arrives at 2R 1 q C logt. Here the central charge 
is c = 1. However just like in Sec. 5.1, subleading corrections require a new and more thorough analysis, as 
other leading irrelevant operators could play a role. This is also an interesting question for further studies. 

6. Conclusion 

We have studied in this paper the universal properties of the emptiness formation probability, mostly focusing 
on a CFT approach but also trying to combine it with other numerical and exact methods. It has been 
known for some time that the EFP can exhibit strikingly different behavior, depending on the model (e.g. 
Ising versus XXZ). We have shown how the Ising universality class can be systematically understood, and 
performed extensive checks of our results in a simple free fermionic realization. It however applies away from 
free fermions, and can also be generalized to any minimal model. For example one could define an analogue 
of the EFP in anyonic chains [95] related to RSOS models [96, 97], and obtain in the same way universal 
terms for this quantity. 

We used extensively the connection with the theory of Toeplitz and Toeplitz+Hankel determinants. 
While the application of such Fisher-Hartwig related results has become increasingly frequent in the physics 
literature, it usually allows to go further than field theoretical calculations. Here we managed to find an 
example where the CFT approach uncovers a £~ 1 \ogt that has yet to be proved. We gave an heuristic 
argument (discontinuities in the derivative of the "regular" part of the symbol) to explain its origin, but 
it would be interesting to perform a more rigorous and systematic study, especially in the more involved 
Toeplitz+Hankel case. In the CFT framework such occurrence requires some fine-tuning: it is a specific 
feature of corner free energies with angle 27r, and shouldn't appear in correlations of local physical observables, 
or in the entanglement entropy. In general both approaches could be pushed much further, and we hope 
that our method can shed some more light on the general structure of subleading corrections to the Fisher- 
Hartwig formula, and its relation to the Renormalization Group in statistical mechanics. As an aside, we 
showed that the I -1 log£ also appears in the full counting statistics for spins in the z basis. 

We also studied the XXZ chain, focusing mostly on the free fermion point. Here the situation is 
somewhat incongruous, as the simplest correlation from the Algebraic Bethe ansatz point of view is actually 
difficult to understand at the field-theoretical level. Our analysis was mostly numerical. We showed that the 
general scaling could be explained in terms of an arctic phenomenon, familiar in the study of dimers, and 
determined the frozen region numerically. We also provided some evidence for the universality of logarithmic 
corrections. A field-theoretical determination of these terms would require a better understanding of the 
fluctuations around the background we identified, and is left as an interesting open question. 

Finally it might also be interesting to investigate similar setups for quantum quenches in real time. 
Indeed the conformal Ising case of Sec. 2 bears some similarity to the Calabrese Cardy argument for global 
quenches [35], where a system is prepared in the ground-state \ipo) of a massive Hamiltonian, but evolves 
at time t > with a critical Hamiltonian. The initial wave function |?/>o) is then seen, in imaginary time, 
as sufficiently close [98] to a conformal boundary state. Interestingly in the XXZ chain it is possible to 
study a slightly different problem where one starts from a strongly magnetized state (see e.g. Ref. [99]), 
which is obviously very far from a conformal boundary state. This type of initial state is exactly the kind 
that generates an arctic phenomenon in imaginary time. It would therefore be interesting to investigate the 
possible influence of this arctic phenomenon in real time quench problems. 
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Appendix A. Numerical determination of the arctic region 



Appendix A. 1. The dimer problem 



Let us consider the dimer problem on a square lattice with open boundary conditions. As is well known, the 
partition function (number of dimer converings) is given by 



Z = ViK = VI detiT| (A.l) 

where K is a signed adjacency matrix. The matrix elements are shown in Fig. Ala. For a given size L x x L y , 
the Kasteleyn matrix can be diagonalized exactly, and its determinant evaluated in closed form. The local 
statistics can be handily accessed using K. For example, the probability to observe a given dimer on the 
link is given, up to a global sign, by 

P(i,j) = Kij x (K-% (A.2) 
K can be inverted for finite L x , L y , we refer to [40] for the exact expressions. In the limit L Xl L y — > oo these 
expressions reduce to integrals that can be computed exactly (see again [40]). We are interested in a system 
where an alternating sequence 1, 0, 1, 0, . . . of dimer occupancies along a segment is imposed somewhere on 
the left of the system, as is shown in Fig. Alb. The number of dimer coverings compatible with such a 
condition is given by 

Z = PfK, (A.3) 

where K is deduced from K by setting to zero the matrix elements corresponding to the removed links in 
Fig. Alb. Inverting K is more difficult, and accessing the puj) should be much harder. There is however a 
simple trick, explained in the next section, that allows to overcome this difficulty. 
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(a) Kasteleyn orientation of the lat 
(i j), Kij — —1 if (i <— j), Kij — otherwise 



ice: Kij — 1 if 



(b) Links removed to ensure the 1,0,1,0 dimer 
occupancies along an horizontal segment. 



Figure Al: (a) Kasteleyn weighting of the lattice, (b) Probability of formation of a "ferroelectric string" of 
length £ = 4 with dimer occupancies 1,0, 1,0. The zigzag lines corresponds to the link we have to remove 
to ensure the ferroelectric configuration. The modified Kasteleyn matrix is obtained by simply setting to 
zero the corresponding matrix elements. The sites part of a removed link, and are represented in red circles. 
There are n — 3£ of them. This will be the set a in Sec. Appendix A.2. 



Appendix A.2. Perturbed Kasteleyn matrices 
K is given by 

K = K + E, (A.4) 

where E is a "perturbation" matrix that has very few non zero matrix elements: — —Kij if the link (i,j) 
has been removed, otherwise. The inverse may be written as 

k~ l = [K + EY 1 (A.5) 

= [^(l + A- 1 ^)]' 1 (A.6) 

= (l + K-^y 1 K' 1 (A.7) 

2G 



Since we know how to calculate K 1 , let us study 

J = K~ 1 K = 1 + K~ 1 E (A.8) 

in more details. Denote by a the set of sites part of a removed link. There are n = SI of them. The other 
sites (the majority) are in b, and there are N — n of them. The matrix J is close to the identity, in the sense 
that it differs from it only for n columns. Indeed, 



(A.9) 



is non zero only if j € a. Now we write J in block form as 

J = 



\ Jba 



1 



bb 



(A.10) 



where a b isanx (N — n) with zero matrix elements, and lbb is the (N — n) x (N — n) identity matrix. 
Using this form, the emptiness formation probability simplifies into 



dctK 
dot if 



de^K^K) = dct J = det J a 



(A.11) 



Such types of formula have already been derived for the monomer correlation function on the square [40] 
and triangular [f 00] lattice. Inverting J, we get 

/ (j^r 1 o ab \ 



J- 1 



(A.12) 



V -4a (4a) lbb J 

Notice ( Jaa) 1 — i-^ 1 ) , an d thus the notation is unambiguous. This is a huge simplification, because 
numerically inverting J actually only requires to invert a small block J aa , the rest following from matrix 
multiplication. Similarly if we write K~ l in block form, we finally get 



k~ 1 = 



■tf (^" 1 ) ab 



The crucial piece of information is encoded in J aa . In the end, the probability of interest reads 



P{ij) = P(i,j) 



(A.13) 



(A.14) 



■^y i-^i \^ aa 

p£a g£a 

e. g. for i € b, j € b. The numerical evaluation of such double sums can be done relatively fast. By that 
we mean that all the probabilities can be obtained in a time of order 0(£ 2 L x L y ), on a L x x L y grid. Let 
us finally mention that inverting J aa is numerically very unstable. Typically, the system sizes shown in this 
paper require close to a hundred digits of precision. 



Appendix A. 3. Exact expressions for the propagators 

The exact expressions are, with i = i x + (i y — 1)L X and j = j x + (j y — 1)L X 
K^ 1 = K- l (i x ,i y ;j x ,j y ) 



4( — 1) 3 ( l +3!c - <• +3y-iy ) 
(L x + l)(L y + l) 



EE 



sin (i x k x ) sin (j x k 

v"yky) sin (j y k y ) 
cos 2 (k x ) + cos 2 (k y ) 



X f{k X , ky) 



(A.15) 
(A.16) 
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with 



and 



cos (k x ) 



(jx - ix) odd, (j y - i y ) even 



f{k x ,k y ) = < coa (k y ) ; (j x - i x ) even, (j y - i y ) odd 



(A.17) 



otherwise 



m x TT 



,m x = 1,2,.. .,1^/2 



L x + 1 
m y TT 



(A.18) 
(A.19) 



In the limit L x , L a — > oo these expressions reduce to integrals that can be computed exactly [40] 
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